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There are two types of asymptotic formulas for the 12 j symbol with one small and 11 large angular 
momenta. We have derived the first type of formula previously in [L. Yu, Phys. Rev. A84 022101 
(2011)]. We will derive the second type in this paper. We find that this second asymptotic formula 
for the 12j symbol is expressed in terms of the vector diagram associated with two 6j symbols, 
namely, the vector diagram of two adjacent tetrahedra sharing a common face. As a result, two 
sets of Ponzano-Regge phases appear in the asymptotic formula. This work contributes another 
asymptotic formula of the Wigner 12j symbol to the re-coupling theory of angular momenta. 

PACS numbers: 03.65. Sq, 02.30.1k, 03.65.Vf 



I. INTRODUCTION 

The Wigner 12 j symbol is described in various text- 
books on angular momentum theory [H-[3]- It has appli- 
cations in atomic physics [!, Q and loop quantum gravity 
Q- In this paper, we derive an asymptotic formula for 
the 12 j symbol in the limit of one small and 11 large 
angular momenta. 

There are two special formulas for the exact 12j sym- 
bol when one of its 12 arguments is zero These spe- 
cial formulas are displayed in Eq. |T]) and Eq. @. In 
an earlier paper Q , we derived an asymptotic formula of 
the 12j symbol where the zero parameter in Eq. ([TJ) is 
replaced by a small parameter, and the other 11 param- 
eters are taken to be large. In this paper, we will derive 
an asymptotic formula for the 12j symbol where the zero 
parameter in Eq. @ is replaced by a small parameter, 
and the other 11 parameters are taken to be large. 

The main theoretical tool we use is a generalization 
of the Born-Oppenheimer approximation, called multi- 
component WKB theory |9h13|]. in which the small an- 
gular momenta are modeled by exact linear algebra, and 
the large angular momenta are modeled by a WKB wave 
function. Each wave function in this model consists of a 
spinor factor and a factor in the form of a scalar WKB 
solution. A gauge-invariant expression for the resulting 
multicomponent WKB wave function is developed in the 
semiclassical analysis of the 9j symbol with small and 
large angular momenta 13]. This gauge-invariant ex- 



pression plays a crucial role in deriving the results in 
ul E3) as wen as the result in this paper. Thus, this 
paper assumes familiarity with it. 

In our earlier paper [9], we find that the first type 
of asymptotic formula for the 12j symbol is based on 



the geometry associated with the 9j symbol on the right 
hand side of Eq. (JlJ. Thus, we expect the second type of 
asymptotic formulas for the 12j symbol with a small an- 
gular momentum to be expressed in terms of geometries 
associated with the two 6j symbols on the right hand 
side of Eq. @. This is in fact the case. The asymptotic 
formula of the 6j symbol in terms of the Ponzano-Regge 
phase UM is well known from the role it plays in Regge 
gravity |15| and topological quantum field theory (l6lll7j. 
Using this asymptotic formula for the 6j symbol, we find 
that the second type of asymptotic formula for the 12 j 
symbol contains two Ponzano-Regge phases. 

We will now give an outline of this paper. In Sec. [TTJ 
we display two special formulas for the exact 12 j symbol, 
and then express the 12j symbol as an inner product be- 
tween two multicomponent wave functions. In Sec. Mil 
we use the procedure outlined in [13[ to find the gauge- 
invariant multicomponent WKB form of these wave func- 
tions. In Sec.|TV] based on the methods developed in [l8j . 
we sketch the semiclassical analysis of these Lagrangian 
manifolds and use the Ponzano Regge formula to confirm 
the result of the analysis. In Sec.fV] we find the spinor in- 
ner products at the intersections of the Lagrangian man- 
ifolds. Putting it all together, we derive an asymptotic 
formula for the 12j symbol in Sec. I VII The last section 
contains comments and discussions. 



II. THE 12j-SYMBOL 

The 12j symbol was first defined by Jahn and Hope Q 
in 1954. The appendix of their paper gives two special 
formulas for the 12 j symbol when one of its 12 parameters 
is zero. In the following, we rewrite Eq. (A9) and Eq. 
(A8) in that appendix using a more convenient labeling. 
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where the square bracket notation [c] stands for [c] — 
2c + 1. In this paper, we will replace the zero parameter 
in Eq. ([2]) by a nonzero small parameter ji = s 2 , and 
take the other 11 parameters to be large compared to s 2 . 

From Eq. (A4) in j|, the 12 j symbol is defined as a 
scalar product 



Jl S 2 312 J346 
33 34 J34 jl35 

ji3 324 is je 
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{ [j'12] [j34] [jl3] [j'24] [7346] [jl35] } 2 

(3) 

where the square bracket notation [•] again denotes [c] = 
2c+l, and \a) and \b) are normalized simultaneous eigen- 
states of lists of operators with certain eigenvalues. We 
will ignore the phase conventions of \a) and \b) for now, 
since we did not use them to derive our formula. In our 
notation, the two states are 



h S 2 -^3 h I5 h J 12 J34 J§46 Jtot 

jl S 2 33 34 35 je jl2 J34 J346 
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(5) 



In the above notation, the large ket lists the operators 
on the top row, and the corresponding quantum numbers 
are listed on the bottom row. The hat is used to distin- 
guish differential operators from their symbols, that is, 
the associated classical functions. 

The states \a) and \b) live in a total Hilbert space of 
six angular momenta Hi <S> H3 <£> H4 <£> ^5 ® 'H.e ® W s , 
where s = s 2 - Each large angular momentum 3 r , r — 
1,3,4,5,6, is represented by a Schwinger Hilbert s pac e 
of two harmonic oscillators, namely, H r = L 2 (R 2 ) [l9| . 
The small angular momentum S is represented by the 



usual 2s + 1 dimensional representation of SU(2), that is, 
H s = C 2s+1 . 

Let us now define the lists of operators in Q and 
([5]). First we look at the operators I r , r — 1,3,4,5,6, 
and J| 4 , J| 46 , J13, J135, which act only on the large 
angular momentum spaces each of which can be 
viewed as a space of wave functions ip(x r i, x r2 ) for two 
harmonic oscillators of unit frequency and mass. Let 



= (x ri i. + iPrn)/V2 and 



jt = 

rfi 



(x 



rfj. 



ip 



77 



)/V2, 



ji = 1,2, be the usual annihilation and creation oper- 
ators. The operators I r and J r i are constructed from 
these differential operators a and as follows, 



I r ^ d^d r , 



(6) 



where i = 1,2,3, and <Ji are the Pauli matrices. The 
quantum numbers j r , r — 1,3,4,5,6 specify the eigen- 
values of both I r and J 2 , to be j r and j r {jr + l)j respec- 
tively. 

The operators J 2 4 , J^g, Ji 3 , Jf 35 that define interme- 
diate coupling of the large angular momenta are defined 
by partial sums of 3 r , 



1 34 



'346 



J 3 + J 4 + J 6 



Jl3— J1+J3: 



Jl35 — Jl+J3 + J5- 
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(8) 



The quantum numbers ji , i — 34, 346, 13, 135 specify 
the eigenvalues of th e op erators Jf to be ji (ji + 1), for i = 
34, 346, 13, 135. See [19( for more detail on the Schwinger 
model. 

Now we turn our attention to the operator S 2 that act 
only on the small angular momentum space C 2s+1 . Let S 
be the vector of dimensionless spin operators represented 



3 



by 2s + 1 dimensional matrices that satisfy the SU(2) Hamiltonians A 
commutation relations the matrix U, 



The spinor r M is the fi th column of 



[Si, Sj] = i £ijk Sk ■ 



(9) 



The Casimir operator, S 2 = s(s + 1), is proportional to 
the identity operator, so its eigenvalue equation is triv- 
ially satisfied. 

The remaining operators J 2 2 , J241 and J tot are non- 
diagonal matrices of differential operators. They are de- 
fined in terms of the operators I r , 3 r i, and Si as follows, 



r^ ] {x,p) = U afi (x,p). 



(16) 



where p = dS(x)/dx. 

The Weyl symbols of the operators I r and j r j, r — 
1,3,4,5,6, are I r — l/2 and J r i, respectively, where 



Ir = 



Jri ~ ~ / ^ &r 'flip r ) (J>v%rv ; (1^) 



(il 2 )a/3 = [Jl + h 2 s(s + l))8 af} + 2Jx • S Q(3 , (10) 

= [Jl + n 2 s(s + i)]s aP + 2J 4 • s aP , (ii) 

(Jtot)a/3 = (Ji + J.3 + J4 + J5 + ^e)S a p + hS a p.(l2) 

These three operators act nontrivially on both the large 
and small angular momentum Hilbert spaces. 



III. 



MULTICOMPONENT WKB 
WAVEFUNCTIONS 



We follow the approach used in [13j to find a gauge- 
invariant form of the multicomponent wave functions 
ipa( x ) — (x,a\a) and ^ b a {x) — (x,a\b). Let us focus 
on t/j^x), since the treatment for ip b is analogous. We 
will drop the index a for now. 

Let £>i, i = 1, . . . , 12 denote the the operators listed in 
the definition of |o) in Eq. (HJ). We seek a unitary opera- 
tor U, such that Di for alii = 1, . . . , 12 are diagonalized 
when conjugated by U. In other words, 



fj, v i 



(13) 



where Aj, i = 1, . . . , 12 is a list of diagonal matrix opera- 
tors. Let 000 be the simultaneous eigenfunction for the 
diagonal entries Xi of the operators Aj, i — 1, . . . , 12. 



Then we obtain a simultaneous eigenfunction 
original list of operators Di from 



00 



of the 



(14) 



Since we are interested in i\) a only to first order in h, all 
we need are the zeroth order Weyl symbol matrix U of U, 
and the first order symbol matrix Aj of A^. The resulting 
asymptotic form of the wave function ip{x) is a product 
of a scalar WKB part Be lS and a spinor part r, that is, 



(15) 



Here the action S(x) and the amplitude B{x) are simul- 
taneous solutions to the Hamilton- Jacobi and the trans- 
port equations, respectively, that are associated with the 



are the 



and where z r ^ — x r ^ + ip Til and z r ^ = x rfl — ip r)Jl 
symbols of d and d^, respectively. The symbols of the 
remaining operators have the same expressions as Eqs. 
©, ©, CEUJ-CES]), but without the hats. 

Among the operators Di, J 2 2 and the vector of the 
three operators J to t are non-diagonal. By looking at (fTU|) , 
the expression for J 2 2 , we see that the zeroth order term 
of the symbol matrix J 2 2 is already proportional to the 
identity matrix, so the spinor r must be an eigenvector 
for the first order term Ji • S. Let t^(Ji) be the eigen- 
vector of the matrix Ji • S with eigenvalue /uJi, that is, 
it satisfies 



(Jl • S) a/3 Tf^ = \x3 x T 



>) 

03 



(18) 



where /i = — s, ... , +s. In order to preserve the diagonal 
symbol matrices J\ through the unitary transformation, 
we must choose the spinor t^' to depend only on the 
direction of Ji. One possible choice of is the north 
standard gauge, (see Appendix A of [Hj]), in which the 
spinor 5 a M is rotated along a great circle from the z-axis 
to the direction of Ji. Explicitly, 



.00 



(Ji) 



) . 



(19) 



where (0i,0i) are the spherical coordinates that spec- 
ify the direction of Ji. Note that this is not the only 
choice, since Eq. ([TBI is invariant under a local U(l) 
gauge transformations. In other words, any other spinor 



r that is related to r by a U(l) gauge trans- 



s ig(Ji 

formation satisfies Eq. (fT8)) . This local gauge freedom is 
parametrized by the vector potential, 



» 



i(rW) 



0Jl 



(20) 



which transforms as A^' = — (g) under a lo- 
cal gauge transformation. Moreover, the gradient of the 
spinor can be expressed in terms of the vector potential, 
(Eq. (A.22) in [111), as follows, 



<9r00 
dJ~i 



Ji x S 



.00 



(21) 
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Once we obtain the complete set of spinors r^', fj, = 
— s, . . . , s, we can construct the zeroth order symbol ma- 
trix U of the unitary transformation U from Eq. (jT6|) . 

Now let us show that all the transformed symbol ma- 
trices of the operators in Eq. (|4]), namely, the Aj, are 
diagonal to first order. Let us write A[D] to denote the 
operator U^DU, and write A[D] for its Weyl symbol. 
First, consider the operators I ri r = 1,3,4,5,6, which 
are proportional to the identity matrix. Using the oper- 
ator identity 

(22) 

we find 



where in the second equality, we have used the reduced 
Lie-Poisson bracket (Eq. (30) in [l9j]) to evaluate the 
Poisson bracket {r, Ji}, and in the third equality, we 
have used Eq. (f2"Tj) for dr/dJi. Notice the term involv- 
ing S in TL, in Eq. (l29l) cancels out the same term in 

A(Ji 2 ) in Eq. (|27|) . leaving us with a diagonal symbol 
matrix 



[A(Ji 2 )]^ = Ji 



(30) 



Taking the square, we obtain 



[A(/ r )]^ = (I r - 1/2)6^ - i«7 * o/1 {Uoau, Ir} , (23) 

where we have used the fact that the symbol of a commu- 
tator is a Poisson bracket. Since U a ^ — is a function 
only of J12, and since the Poisson brackets {Ji, I r } = 
vanish for all r — 1, 3, 4, 5, 6, the second term in Eq. (|2"3")l 
vanishes. We have 



[A(I r )W = & - 1/2) <W ■ 



(24) 



Similarly, because {Ji, J| 4 } = and {J i, J^g} = 0, we 
find 



[A (■'34)] A** 1 — ^34 <W ' ^346)]^ — ^346 <W • (25) 



Now we find the symbol matrices A(J 12 ) for the vector 
of operators J 12, where 

[A(Ji 2 )W = U^ixSc^Upu + hUl fl S a pUp„ . (26) 

After converting the above operator equation to Weyl 
symbols, we find 

[A(Ji 2 )W (27) 

= Jl^i/ - iHUapiUap, Ji} + hUapSapUpp 

= - »RtW*{tW, Ji} + hr^*S a pT^ . 

Let us denote the second term above by T^„, and use 
?T]), the orthogonality of t, 



O)* _(«0 _ x 

'a <a ~ u fw 1 



to get 



(28) 



(29) 




[A(J? 2 )W = ( J i +/ifi) 2 <W 



(31) 



Finally, let us look at the last three remaining opera- 
tors Jtot in Eq. (fT"2"|) . Since each of the the symbols 
J r for r — 3,4,5,6 defined in Eq. (|17[) Poisson com- 
mutes with J 1; that is, {Ji, J r } = 0, we find A(J r ) = 
J r — ihUl{Ua(3i),J r } = 3 r , for r = 3,4,5,6. Using 
A(Ji 2 ) from Eq. (|3"D]). we obtain 



[A(Jtot)]/i^ 

M 1+ £ 



SA^ x Ji + (J 3 



J 6 ) 



(32) 

8uv ■ 



Therefore all A*, i = 1, . . . , 12 are diagonal. 

Not counting the trivial eigenvalue equation for S 2 , we 
have 11 Hamilton- Jacobi equations associated with the 
Aj for each polarization fi in the 20 dimensional phase 
space C 10 . It turns out that not all of them are func- 
tionally independent. In particular, the Hamilton- Jacobi 
equations A(j 2 ) = J 2 h = (ji + l/2)K and A(Jf 2 ) = 
(Ji + fih) 2 = (ji2 + l/2) 2 h 2 are functionally dependent. 
For them to be consistent, we must pick out the polar- 
ization [i — ji2 — ji . This reduces the number of inde- 
pendent Hamilton- Jacobi equations for S(x) from 11 to 
10, half of the dimension of the phase space C 10 . These 
ten equations define the Lagrangian manifold associated 
with the action S(x). 

Now let us restore the index a. We express the multi- 
component wave function tp® (x) in the form of Eq. (|15[) , 



r a {x) = B a {x)e lS ^/ h T a a {x,p). 



(33) 



Here the action S a (x) is the solution to the ten Hamilton- 
Jacobi equations associated with the /i th entries Xf of 10 
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of the symbol matrices A" , given by 

h = (j 1 + i/2)h, 
h = ( j3 + i/2)h, 
h = {h + i/2)h, 

h = Cj 8 + 1/2)», 
h = Cj6 + l/2)», 
= 0*34 + l/2) 2 ^ 2 , 

\/2) 2 h 2 , 



(34) 



346 
r(o) 



(.7346 
Jl 



1 + ^ 
Jl 



h Ai x Ji + (J 3 + J 4 + J 5 + J 6 



and r a = r^ 1 ) with it = 



Note that all the Hamil- 



tonians in Eq. (|3"4"|) except the last three, J^, preserve 
the vector value of Ji and J 5 along their Hamiltonian 
flows. 

We carry out an analogous analysis for The 
result is 



1 p b a (x)=B b (x)e iS ^ h r b (x,p). 



(35) 



where Sb(x) is the solution to the following 10 Hamilton- 
Jacobi equations: 



where the action S^ s (x) is the integral of p dx starting 
at a point zq, which is the lift of a reference point xq in 
the Lagrangian manifold C G J S . The Lagrangian manifold 
C Q J S is defined by the following equations: 



(40) 



h 


= (h 4 


i/2)h, 


h 


= (J3 4 


1/2)71, 


h 


= (j4 4 


1/2)71, 


h 


= (35 4 


i/2)h, 


h 


- (36 4 


1/2)71, 


I 2 

J 34 


= (j34 - 


f i/2) 2 h 2 , 


J 346 


= (.7346 


+ i/2) 2 h 2 , 


Jtot 


= Ji4 


J 3 + J 4 4 J 5 4 J 6 = 



The rotation matrix U a (x) that appears in Eq. (|39l) is 
determined by the SO (3) rotation that transforms the 
shape configuration of Ji and J5 at the reference point 
zq = (xo,p(xo)) on L G J S to the shape configuration of Ji 
and J5 at the point z — (x,p(x)) on C 6 J S . Here Ji and 
J5 are functions of z and are defined in Eq. (|17[) . 

Similarly, the multicomponent wave function for the 
state I b) has the following form, 



h = (ji - 
^3 = (h ■ 
h = {34 - 

h = (35 - 
h = (k - 
J\z = (m 

•A 2 35 = (31 
r(b) _ 



35 



l/2)H, 

i/2)h, 
1/2)71, 
i/2)h, 
i/2)h, 
-i/2) 2 h 2 , 
■f i/2) 2 h 2 
vh 

+ x 



(36) 



H A 4 x J 4 + (Jj + J 3 + J 5 + J 6 



Here the spinor r b = Tj satisfies 



(J 4 



(37) 



where 1/ = j 2 4 - J4- 

The vector potential A 4 is defined by 



(38) 



Again, note that all the Hamiltonians except the last 
three, J^t, preserve the value of J4 and J 6 along their 
Hamiltonian flows. 

We follow the procedure described by the analysis pre- 
ceding Eq. (69) in [13] to transform the wave functions 
into their gauge-invariant form. The result is a gauge- 
invariant representation of the wave function, 



r(x) = B a {x)e< ia ^ h [U a {x)T a {x Q )\ 



(39) 



4> b (x) = B b (x)e< 33 ^' h [U b (x)T b (x )] 



(41) 



where the action S b ° s (x) is the integral of p dx starting at 
a point that is the lift of xq onto the Lagrangian manifold 



-6js 



The Lagrangian manifold C b s is defined by the 



following equations: 



h = 

h = 

h = 

h = 

h = 
j 2 - 



J 



135 
Jtot 



(Ji4 
(33 1 
(34^ 
0*5 + 
(je + 

(ii3 

(.7135 

Ji 4 J 3 



1/2)71, 
l/2)h, 

i/2)h, 
i/2)h, 

1/2)*, 

i/2) 2 h 2 , 
f \/2) 2 h 2 . 



(42) 



0. 



The rotation matrix U b (x) that appears in Eq. (|41l) is 
determined by the SO (3) rotation that transform the 
shape configuration of J4 and Jg at the reference point 
z o = (%o,p(xo)) on C b JS to the shape configuration of J4 
and Jg at the point z = (x,p(x)) on £ b Js . 

Taking the inner product of the wave functions, and 
treating the spinors as part of the slowly varying ampli- 
tudes, we find 

(b\a) = e 4K J2 n " exp{»[5^ s (z fe ) - S^ s (z k ) - ^/2}/h) 

k 

(uS k T b (z )y(U a a k T*(z Q )). (43) 
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In the above formula, the sum is over the components 
of the intersection set M. k between the two Lagrangian 
manifolds C^ s and C^ s . The point z k is any point in 
the Ath component. The amplitude ilk and the Maslov 
index fi k are the results of doing the stationary phase 
approximation of the inner product without the spinors. 
Each rotation matrix U® h is determined by a path <-y a ( ofc ) 
that goes from zq to Zk along C^ s , and U® k is similarly 
defined. The formula (|43|) is independent of the choice 
of z k , because any other choice z' k will multiply both 
and U° j by the same additional rotation matrix which 
cancels out in the product (U® yU® . 

IV. THE LAGRANGIAN MANIFOLDS 

We now analyze the Lagrangian manifolds C® JS and 
C b JS , defined by the Hamilton- Jacobi equations Eq. (|4TJ)) 
and Eq. (|42l) . respectively. We focus on C 6 J S first, since 
the treatment for C, 3S is analogous. Let it : $5^ — > A$j 
denote the projection of the large phase space $5^ = 
(C 2 ) 5 onto the angular momentum space A$j = (R 3 ) 5 , 
through the functions J r ;, r = 1,3,4,5,6. The first six 
equations, = j r + 1/2, r = 1,3,4,5,6 fix the lengths 
of the five vectors |J r | = J r , r = 1,3,4,5,6. The three 
equations for the total angular momentum, 

Jtot = Ji + J3 + J 4 + J5 + Je = , (44) 

constrains the five vectors J; , i = 1 , . . . , 6 to form a close 
polygon. The remaining two equations 




FIG. 1. The configuration of a point on C^ s , projected onto 
the angular momentum space Asj, and viewed in a single copy 
ofR 3 . 



itself is an Abelian group of eight elements, 
(Z2) 3 = {e, x, y, z, xy, xz, yz, xyz}. Thus the mani- 
fold C 6 J S is topological^/ U(l) 7 x SU(2)/(Z 2 ) 3 . The 
analysis for C b JS is the same. 

Now it is easy to find the invariant measure on C^ s and 
£ 6 b js . It is dipi A dips A dip 4 A dips A dtp 6 A d6> 34 A d9 34e A du, 
where du is the Haar measure on SU(2). The volumes 
Va of C & J S and Vb of C Q J S with respect to this measure 
are 



J 3 2 4 = (J34 + l/2) 2 ^ 2 , (45) 
4 i6 = 4s = 0'846 + l/2) 2 ^ 2 , (46) 

put the vectors J3 , J4 into a 3-4-34 triangle, and put the 
vectors Ji,Js into a 1-5-346 triangle. Thus, the vec- 
tors form a butterfly shape, illustrated in Fig. [TJ This 
shape has two wings (J3, J4, J34) and (Ji, J5, J346) that 
are free to rotate about the J34 and J 3 46 edges, respec- 
tively. Moreover, the Hamilton- Jacobi equations are also 
invariant under an overall rotation of the vectors. Thus 
the projection of C e J s onto the angular momentum space 
is diffeomorphic to U(l) 2 x 0(3). 

The orbit of the group U(l) 5 generated by I r , 
r = 1, 3, 4, 5, 6 is a 5-torus. Thus C®J S is a 5-torus bundle 
over a sub-manifold described by the butterfly configura- 
tion in Fig.[TJ Altogether there is a U(l) 7 x SU(2) action 
on C^ s . If we denote coordinates on U(l) 7 x SU(2) 
by {ipi,ip2, ipz,ip4, ip& ,#34, #346, u), where u g SU(2) 
and where the five angles are the 47r-periodic evolution 
variables corresponding to {I\ , I3 , J4, 1$ , Jg , J34, J346) , 
respectively, then the isotropy subgroup is generated by 
three elements, say x = (27r, 27r, 2tt, 2tt, 2tt, 0, 0, — 1), 
y = (27r,0,0,27r,27r,27r,0, -1), and z = 
(2tt, 0, 0, 27T, 0, 0, 2-7T, — 1). The isotropy subgroup 



V A = V B = - (4tt) 7 x 16tt 2 = 2 15 tt 9 , (47) 
8 

where the 1/8 factor compensates for the 8-element 
isotropy subgroup. 

We now examine the intersections of C^ s and C G b ^ s in 
detail. Because the two lists of Hamilton- Jacobi equa- 
tions Eq. (|40|) and Eq. (|42|) share the common equations 
I r = j r + 1/2, r = 1,2,3,4,6, the intersection in the 
large phase space $5^ is a 5-torus fiber bundle over the 
intersection of the projections in the angular momentum 
space A 5 j. The intersections of the projections in A 5 j 
require the five vectors J r , r — 1, 3, 4, 5, 6, to satisfy 

I Jr I — Jr 5 ^ Jr — , 

r 

|J3+J 4 | = J 3 4, |J 3 + J4 + Je| = ^346, (48) 

|Jl+J 3 | = Jl3, |Jl+J 3 +Js| = Jl35- 

These conditions imply that the six edges 
^3,^135,^346,^4,^34,^6 form a tetrahedron, and 
the six edges J3, J135, J346, Ji, J13, J5 form another 
tetrahedron. The two tetrahedra share the common 
face (J3, J135, ^346), as illustrated in Fig. [5] We can 
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use the procedure explained in the appendix of [20| 
to construct a tetrahedron with the six edge lengths 
J?>, J135, ^346 1 J4, ^34 j J&- This procedure gives us the 
vectors J3, J135, J346, Ji, J5, Ji3- Then we use the 
following three conditions 



J4 • J4 
J 



J A 



I ' J 3 - 2^34 _ J i - J3) 



(49) 



J4 ' Jl35 — ^(^4 + ^135 _ J&) 



to solve for the vector J4. In general, there are two solu- 
tions for J4 , corresponding to the two orientations of the 
second tetrahedron. See Fig. [2] and Fig. [3l Once we have 
J4, we then use 



J34— J3+J4 
J 6 = J346 — J34 



(50) 
(51) 



to construct the second tetrahedron with the six edges 

Jzi ^135) ^346) J4, ^34, Je- 

The two solutions to Eq. (|4T)1) give rise two vector con- 
figurations that are not related by an 0(3) symmetry. 
This means that there are four vector configurations sat- 
isfying Eqs. (|4"5)l that are not related by an SO (3) sym- 
metry. The intersections in $5^ are the lifts of the in- 
tersections in A 5 j. Therefore, the intersection of C^ s 
consists of four disconnected subsets, where each subset 
is a 5-torus bundle over SO(3). Let us denote the two 
sets corresponding to the configuration in which the two 
tetrahedra are on opposite sides of the ( J3, J135, J346) 
triangle and its mirror image by Jn,/i2, and denote the 
configuration in which the two tetrahedra are on the same 
side of the triangle (J3, J135, J346) and its mirror image 
by I21 , 122 ■ The vector configuration for a typical point in 
In is illustrated in Fig. [5J and the vector configuration 
for a typical point in J21 is illustrated in Fig. [3] Each 
intersection set is an orbit of the group U(l) 5 x SU(2), 
where U(l) 5 represent the phases of the five spinors and 
SU(2) is the diagonal action generated by Jtot- 

The isotropy subgroup of this group action is Z 2 , gen- 
erated by the element (2tt, 2-7T, 27r, 27t, 2tt, —1), in coordi- 
nates (ipi,ip3,'ip4,ip5,ip6 7 u) for the group U(l) 5 x SU(2), 
where u G SU(2). The volume of the intersection man- 
ifold In, I12, I21, or I%2, with respect to the measure 
dipi A dip3 A dip 4 A dip5 A dips A du, is 



1 



Vi = ^(47r) b x 16^ = 2 13 tt 7 



(52) 



where the 1/2 factor compensates for the two element 
isotropy subgroup. 

The amplitude determinant is given in terms of a deter- 
minant of Poisson brackets among distrinct Hamiltonians 
between the two lists of Hamilton- Jacobi equations in Eq. 
(HQ)) and Eq. (|4*21 . In this case, those are (J34, J346) from 




FIG. 2. The configuration of a point on the intersection set 
In, projected onto the angular momentum space As-,-. 




FIG. 3. The configuration of a point on the intersection set 
I21, projected onto the angular momentum space As-,-. 



H). Thus the deter- 



Eq. flU and (J 13 , J 135 ) from Eq. 
minant of Poisson brackets is 

{J34, J13} {J34, J135} 
{J346, J13} {<^346, ^135} 



1 



•/34^346<^13>^135 
1 

>/34>/346^13^135 



^341 
Vl35 



V 



:S l(i 



V, 



3(46)(15) 



1^135^346 I 



(53) 



where, in the last equality, we have used V-^is) (46) = 
J 3 ' (J46 x J 15) = Oj since the edges ( J3, J46, J15) form a 
triangle. Here is six times the volume of the tetra- 
hedron generated from Ji , Jj , J/- , and is given by 
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Vi 



ijk 



J 8 • (Jj X J k ) 



(54) 



The amplitude fi/. in Eq. (|4"3")l can be inferred from Eq. 
(10) in [21j. In the present case, each fl k has the same 
expression il. It is 



n 



[2m)Vl V ^34^346^13^135 



VTTO ^1^135^3461 
(27Ti)2 13 7r 7 %/J 3 4J346^13>/l35 



2 15 7T 9 



VI ^135 ^346 I 



*V ^34</346^13^: 



135 



2^1^135^346 1 



(55) 



We now outline the calculation of the relative phase 
between the exponents S a (z\2) — 6*6(212) and S a (zu) — 
Sb{zu), which can be written as an action integral 



= (S a (z 12 )-S b (z 12 ))-(S a (z u )-S b (z n )) = * pdx 



(56) 

around a closed loop that goes from z\\ to z\i along C® JS 
and then back along £ b Js . 

We shall construct the closed loop giving the relative 
phase by following the Hamiltonian flows of various 
observables. This loop consists of four paths, and it is 
illustrated in the large phase space <fr 5 j in Fig. 2) The 
loop projects onto a loop in the angular momentum space 
A$j , which is illustrated in Fig. [5] We take the starting 
point p £ In of Fig. @]to lie in the 5-torus fiber above 
a solution of Eq. (|48l) . The projection of p in A$j is 
illustrated in part (a) of Fig. [5j 




FIG. 4. The loop from a point p 6 In to q £ J12 along C e J" 
and then to q' £ I12 along I12, and then to p' 6 Jn along 

and then to p along 7n. 



jC^' 3 , and finally back to 



First we follow the J§ 4 -flow and then the Jj^g-flow to 
trace out a path that takes us along C^ s from a point p in 
In to a point q in Ii 2 . Let the angles of rotations be 2^34 
and 20346, respectively, where </>i2 is the angle between 
the triangles 3-4-34 and 34-6-346, and 034g is the angle 
between the triangles 34-6-346 and 1-5-346. These rota- 
tions effectively reflect all five vectors J r , r = 1, 3, 4, 5, 6 





\A4 






Jg 




(a) 












r) \ ^^^^^ ^ 346 










/■h 







h \ \ | 
\ ^ - - * 






\-h / 




1 /./4 


(b) 






(c) 



(d) 



FIG. 5. The loop from Fig. [4] projected onto a loop 
(a)-Kb)-Kc)->-(d)-Ka) in A 5j . 



across the triangle 1-5-346, taking us from part (a) to 
part (b) of Fig. |5J 

Next, we follow the Hamiltonian flow generated by 
— J346 ■ Jtot along /12, which generates an overall rota- 
tion of all the vectors around — J346- Let the angle of 
rotation be 2034g defined above. This brings the triangle 
34-6-346 back to its original position. However, the tri- 
angle 1-5-346 is now rotated to the other side of triangle 
34-6-346, as illustrated in part (c) of Fig. [5] Effectively, 
the actions on all five vectors 3 r , r = 1,3,4,5,6 from 
part (a) to part (c) of Fig. [5] has been to reflect them 
across the 34-6-346 triangle. 

To go back to a point p' in In, we follow the Jf 3 - 
flow and J^ 35 -flow along CJ S . Let the angle of rotations 
be 2^13 and 2^135, respectively, where 0i3 is the angle 
between the triangle 1-3-13 and the triangle 13-5-135, 
and 0i35 is the angle between the triangle 13-5-135 and 
the triangle 4-6-135. These rotations effectively reflect 
all the vectors across the 4-6-135 triangle, taking us from 
part (c) to part (d) of Fig. O Thus we arrive at a point 
p' € In- We now use the fact that the product of two 
reflections is a rotation about the intersection of the two 
reflection planes. We note that the vector Jg is stationary 
under the reflection across the 34-6-346 plane, as well as 
across the 4-6-135 plane. Thus, the final rotation that 
brings all the vectors back to their original positions is 
generated by — jg ■ Jtot along In. It is an overall rotation 
of all the vectors about — j 6 by an angle 20 6 , where 0g 
is the angle between the 34-6-346 triangle and the 4-6- 
135 triangle. This rotation takes us from part (d) back 
to part (a) of Fig. [5] We denote the final point by p" 
in the large phase space. The points p and p" have the 
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same projection in the angular momentum space A§j. 
Thus the two points p and p" differ only by the phases of 
the five spinors, which can be restored by following the 
Hamiltonian flows of (Ji, ^3, 14, 15, I§). This constitutes 
the last path from p" to p. 

To summarize the rotational history in the angular mo- 
mentum space, we have applied the rotations 



-R(-j6, 20 6 )-Rl35(jl35> 2 ^135)-Rl3(ji 3l 20i 3 )-R(-j346, 20346) 



-R346(j346, 20346^34 (j34, 2034) , 



(57) 



To undo this pure phase, we follow the Hamiltonian flow 
of I§ by an angle —206. 

Similarly, we can find the rotations acting on Ji , J3, J5, 
and proceed to calculate the action integral as in [18] . 
Instead, we will take a shortcut and use the fact that 
the two Lagrangian manifolds C Q J S and C b ls describe the 
WKB wave functions associated with the product of two 
6j symbols on the right hand side of Eq. @. The asymp- 
totic limit of a product of two 6j symbol can be easily 
derived from the Ponzano-Regge formula for a single 6j 
symbol. 



where i? 3 4 acts only on J 3 and J4, -R346 acts only on J 3 , 
J4, and Jg, -R13 acts only on Ji and J3, -R135 acts only 
on Ji, J 3 , and J 5 , and i?(-j 3 46, 20346), -R(-J6,20 6 ) acts 
on all five vectors. The corresponding SU(2) rotations, 
with the same axes and angles, take us from point p in 
Fig. 0] to another point p" along the sequence p — > q — » 
q' -> p' -> p". 

To compute the final five phases required to close the 
loop, we use the Hamilton- Rodrigues formula |22[, in the 
same way as Eq. (46) in [18]. Let us start with vector J4. 
The action of the rotations on this vector can be written 



J346 J3 J135 
h k hi 
1 

27rVl^l35^34 6 | 
1 



J346 J3 J135 

J13 is h 



(60) 



os (s 2 + I) cos (flf x + \) 
cos (Si + S 2 + |j + cos (Si - S 2 ) 



47r V / |V 135 V r 346| 
where the Ponzano-Regge phases Si and S 2 are 



Si = J 4 04 + ^6^6 + ^34034 + J 3$$ + ^135^135 + 0346 0346 



(61) 



206)-R(-j346, 20346) 
#3460346) 2034 6 )i? 3 4 (j 34 , 20 34 )J 4 
= ^(-J6, 20 6 )i?34(j34, 2034)J4 
= i?(j 4 ,20 4 )J4 

= J 4 • (58) 

where we have used the Hamilton-Rodrigues formula in 
the second equality. Thus, we find that the product of the 
rotations acting on J4 is i?(j4, 204), where 04 is the angle 
between the triangle 3-4-34 and the triangle 4-6-135. We 
can lift the rotations up to SU(2) with the same axis and 
angle. Its action on the spinor at p is a pure phase. To 
undo this pure phase, we follow the Hamiltonian flow of 
I4 by an angle —204, modulo 2tt. 

For the vector J6, the rotations acting on it is simple. 



R(— J6; 206)i?(— j346, 2034 6 )i?34 6 (j 3 4 6 , 20 3 4 6 )J 6 

= i?(-j 6 ,20 6 )J 6 

= J 6 (59) 

We can lift the rotations up to SU(2) with the same axis 
and angle. Its action on the spinor at p is a pure phase. 



S2 = Jl01 + Js05 + JlS^W + ^303 + ^135"0135 + ^346 0346 

(62) 

Here -04, 06j 034j 0135, 0346 are the exterior dihe- 
dral angles of the tetrahedron formed by the six edges 

Ji, J e , J34, J3, J135, 0346, and 0i, 5 , 0i 5 , 0^ 2) , 0^, ip$ 6 
are exterior dihedral angles of the tetrahedron formed by 
the six edges Ji, J5, J15, J3, J135, J346- These two tetra- 
hedra are illustrated in Fig. [2l 

The action integral J pdx along the loop p — » q — » 
q' — > p' — > p" — > p in Fig. [4] can then be read off from Eq. 
It is given by 

= 2(Si + S 2 ) . (63) 

The action integral along a similar closed loop from I 2 i 
to I 22 and back to J21 is 

S { ^ = 2(Si - S 2 ) . (64) 

The Maslov indices fii = — 2 and /i 2 = can also be read 
off from Eq. (|6"0|) . Putting the amplitudes fl from Eq. 
([53)1 and the relative actions S^ and S^ 2 ' and Maslov 
indices /ii and fx 2 into Eq. (|4"3")l . we find 
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(%) = ^#S^{e- [(r 6 (z 11 ))t(r«(, 11 )) + e^ (I) +^/^ (llfV( ftl )) t feM*i))" 

■«r A/ 1 V135 V346 1 L V ' X ' . 

+e** [(r fc (z 21 )) t (r a (z 21 )) + e l(S<2,)/R ^ r b (z 21 )f (ttfV(»i))l } 



(65) 



Here we have factored out two arbitrary phases e lKl and 
e iK 2 f or the two pairs of stationary phase contributions. 

The rotation matrices Ud , i = 1, 2 are determined by the 
paths from zn to Zi2 along L e J s . Similarly the rotation 

(i) 

matrices Uj. , , i = 1, 2, are determined by the paths from 
zn to z i2 along C 6 b JS . 

V. THE SPINOR PRODUCTS 

We now calculate the spinor products in Eq. ff65|) . We 
choose the vector configurations associated with z\\ to 
correspond to a particular orientation of the vectors. We 
put Ji along the z-axis, and put J5 inside the xz-plane, 
as illustrated in Fig. [51 Let the inclination and azimuth 
angles (9, <fi) denote the direction of the vector J4. From 
Fig. [6l we see that 4> is the angle between the (Ji, J5) 
plane and the (Ji,J4) plane. We denote this angle by 
4> = 4>i ■ The inclination angle 9 is the angle between the 
vectors Ji and J4. 

Z 

1 1 




X 



FIG. 6. The vector configuration at the point zn in In. 



T a a {zn) = <5 QM . (66) 

For the spinor T b (zn), we choose it to be an eigenvector 
of J4 • S in the north standard gauge, that is, 

T b a (z n ) = e« a -^ dl a (9) . (67) 
Taking the spinor inner product, we obtain 

(r b (z 11 ))t(r a (z 11 )) = e ^-»)* d ^ {e) . ( 68 ) 

To evaluate the other spinor product at Z12, we need 
to find the rotation matrices and , which are 
generated from paths j a and 7^ from zu to Z12 along 
C^ s and C b ls , respectively. 

We choose the path j a to be the path from p to q 
generated by the J| 4 -flow and the J§ 46 -flow, which are 
illustrated in Fig. [4] in the large phase space, and in part 
(a) of Fig. [S]m the angular momentum space. This path 
contains no flow generated by the total angular momen- 
tum, so 

C/W=l. (69) 

we choose the path 7^ to be the inverse of the path from q 
back to p along C G J S in Fig. 0] which contains the overall 
rotations 

U ( b 1] = tf(j 3 46, 20346 Mk 2^ 6 ) ■ (70) 

Because only overall rotations can move the vectors J4 
and J6 along the flows on C b s , we can determine this 
rotation by looking at its effect on J4 and Jq. The effect 
of the rotation on J4 and 3q is to reflect them across the 
1-5-346 triangle. In the particular frame that we chose, 
the rotation U b effectively moves J4 to its mirror image 
J 4 across the 1-5-346 triangle in the xz-plane, which has 
the direction given by (—(pi, 6). Thus UbT b (zu) is an 
eigenvector of J 4 • S, and is up to a phase equal to the 
eigenvector of J 4 • S in the north standard gauge. Thus, 
we have 



The gauge choices for the spinors at the reference point 
zn are arbitrary, and they only contribute a phase that 
can be absorbed into e iKl . To be concrete, since Ji points 
in the z direction, we choose the spinor r°(zn) to be the 
fi th standard eigenvector for S z , that is, 



[Ul 1] r b (z u )] Q = e ivH * e -^ a -^ d s va {9) , (71) 

where H4 is a holonomy phase factor equal to the area of 
a spherical triangle on a unit sphere. See Fig. [7] There- 
fore, the spinor product at the intersection I12 is 
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(U { b 1} r b (z 11 ))t(C/ Q r a (z 11 )) = e ivH * d a u J9) . 



(72) 




FIG. 7. The phase difference between two gauge choices can 
be expressed as an area around a closed loop on the unit 
sphere. 



Let us denote the first term in Eq. (|65|) by T\. Substi- 
tuting the spinor inner products Eq. ([68| and Eq. (|72|) 
into Eq. (|65[) , we find that T\ is given by 



- r ini V^34J346./l3./l35 , s f(j s 

1 i=e ff= = I 



cos 



7T a/ I ^135^346 1 

2 



(73) 



- + fi4>i + v l — - 0i 



Using a different choice of the reference point in Fig. 
(|SJ) and a different set of paths, we can derive an alter- 
native expression for the inner product, and eliminate 
the term H^. Let us choose a new reference point zu 
to correspond to an orientation in which J4 is along the 
z-axis, and Jg lies in the x-z plane. Through essentially 
the same arguments as above, we find 



r wt! V ^34^346^13^135 „ , Q s 



COS 



7T\/ 1 ^135^346 1 



(74) 



— + 2 + "Pi + V 



Hi 

2 



Here Hi is another holonomy for the Ji vector, and the 
angle 4>4 is the angle between the (J4, Jg) plane and (J4, 
Ji) plane. Because the quantities 4>i, 04, Hi,H^ de- 
pend only on the geometry of the vector configuration, 
and are independent of /i and v, we conclude that the ar- 
gument in the cosine must be linear in fi and v. Equating 
the two arguments of the cosine in Eq. (1731 and in Eq. 
(1741 . we find that this linear term is (/x0i + V(f>4). We 
find 




FIG. 8. A different choice of vector configuration at zn in 
hi- 



4l V J34 </346>/l3>/l35 
7T\/|Vi35V346| 



(75) 



,(1) ,(1) 7T 

— + M ] + vA ] + 2 



Here we have put back the indices (1) to indicate that 
we are using the vector configuration in which the two 
tetrahedra are on opposite sides of the 3-346-135 triangle. 
Through an analogous calculation, we find 



T 2 = e l 



V >/34>/346^13^135 



COS 



"V 1^135^346 1 

— - m + 



(76) 



A 2 ) 



Here the indices (2) indicate that we are using the vector 
configuration in which the two tetrahedra are on same 
side of the 3-346-135 triangle. 



VI. AN ASYMPTOTIC FORMULA FOR THE 
12j-SYMBOL 

From the definition Eq. ([3]), we see that the factor 
([J34][j346][ji3][ji35]) 1/2 in the denominator of Eq. © 
partially cancels out the factor (J34 J346-/i3>/i35) 1 ^ 2 from 
Ti and T2 in Eq. ([75)1 and Eq. ([TBI , respectively, leaving 
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a constant factor of 1/4. Because the 12j symbol is a real 
number, the relative phase between e IKl and e iK2 must 
be ±1. Through numerical experimentation, we found it 
to be (— l) 2s2 . We use the limiting case of 32 = s = 
from Eq. ([2]) to determine the overall phase convention. 



This determines most of the overall phase. The rest can 
be fixed through numerical experimentation. Putting the 
pieces together, we obtain an asymptotic formula for the 
12j symbol with one small quantum number: 



Jl S J12 Jl25 
33 34 J34 3 135 
3 13 J24 35 36 



4^^/1^135^3461 



{d°{e {1) ) 



cos 



Si +S 2 



(77) 



-(-1) 2 X^ (2 



S\ — S2 



A 2 ) 



A 2 ) 



10- 8 




FIG. 9. Comparison of the exact 12j symbol (vertical sticks 
and dots) and the asymptotic formula (|77[) in the classically 
allowed region away from the caustics, for the values of j's 
shown in Eq. (fg2j) . 



ht 10 




FIG. 10. Comparison of the exact 12j symbol (vertical sticks 
and dots) and the asymptotic formula Q77p in the classically 
allowed region away from the caustics, for the values of j's 
shown in Eq. (p3)l . 



Here, the indices on the d-matrix are given by /i = 
3 12— ji and v = 324—34- They are of the same order as the 
small parameter s. The two Ponzano-Regge phases Si 
and S2 are defined in Eq. (|pTj) and Eq. (lrJ2l , respectively. 
The quantity Vijk is defined by 



(78) 



The angle 9 is the angle between the vectors J12 and 
J 13. The angles (f>i,<f>4,6 are given by the following equa- 
tions 



7T — cos 



(Jl X J 4 ) ■ (Ji X J 5 ) 

|Jl X J 4 | |Ji X J 5 | 
(J 4 X Ji) • (J 4 X J 6 ) 



IJ4 X JJ |J 4 X J 6 



= cos 



Ji J 4 

J1J4 



(79) 
(80) 
(81) 



These angles are calculated using the vector configu- 
ration at a point in /2I) which is illustrated in Fig. [3] 

We illustrate the accuracy of the approximation Eq. 
(|77|l by plotting it against the exact 12j symbol in the 
classically allowed region for the following values of the 
j's: 



Ji s 2 312 J125 

33 34 334 3135 
3X3 324 35 J6 



35 1 34 39 

36 28 38 31 
27 29 j 5 36 



(82) 



The result is shown in Fig. [9] 

Since the asymptotic formula (|77|) should become more 
accurate as the values of the j's get larger, we plot the 
formula against the exact 12j symbol for another exam- 
ple, 
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3i *2 312 3125 ) { 177/2 5/2 88 89 ) 
33 34 334 Ji35 \ = I 181/2 141/2 87 77 L 
3i3 324 35 k J [ 75 73 j 5 91 J 

(83) 

in the classically allowed region away from the caustic 
in Fig. [lOj These values of the j's are roughly 2.5 times 
those in Eq. (|82j) . Again, the agreement has improved. 

VII. CONCLUSIONS 

In this paper, we have derived the second asymptotic 
formula for the Wigner 12 j symbol with one small angu- 
lar momentum. Eq. ((77)) here and Eq. (80) in [9] together 
cover all the different placements of the small angular 
momentum among the 12 parameters of the 12j symbol. 
Although the two asymptotic formulas of the 12 j symbol 
are similar, in that they are expressed in terms of the 
asymptotic phases of lower 3nj-symbols, Eq. (1771) in this 
paper is in many ways simpler. The relationship between 
the 6j symbol and the geometry of a tetrahedron is well 



known. The construction of the vectors of the tetrahe- 
dra is simpler than that for the vector diagram of a 9j 
symbol in Q. In any case, the two formulas are valid 
even when we take s to be large, as long as the other 11 
angular momenta j are much larger relative to s, that is, 
as long as 1 << ,s << j. 

Currently, the asymptotic formula for the 12 j symbol 
when all the angular momenta are large is still unknown. 
Such a formula must reduce to Eq. (|77| here or Eq. (80) 
in [9j in the limit 1 << s << j. Therefore, the work 
in this paper may eventually help us find the asymptotic 
formula of the 12 j symbol when all j are large. 

After the completion of this manuscript, I was in- 
formed by the authors of [23[ that they have found an 
independent derivation for the asymptotic formula of the 
9j symbol in (T3 |. 
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